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The optical absorption spectrum of the three most stable isomers of the Agn system was calculated using 
the time-dependent density functional theory, with the generalized gradient approximation for the exchange and 
correlation potential, and a relativistic pseudopotential parametrization for the modelling of the ion-electron 
interaction. The computational scheme is based on a real space code, where the photoabsorption spectrum is 
calculated by using the formalism developed by Casida. The significantly different spectra of the three isomers 
permit the identification of the ground-state configuration predominantly present in the laboratory beams in 
base to a comparison between the calculated photoabsorption spectrum of the most stable configuration of Agn 
and the measured spectra of medium-size silver clusters trapped in noble gas Ar and Ne matrices at different 
temperatures. This assignment is confirmed by the fact that this isomer has the lowest calculated energy. 

PACS numbers: 77.22.-f,61.46.Bc,31.15.Ew 

Keywords: silver, cluster, photoabsorption spectrum, density functional theory 



I. INTRODUCTION 

The study of transition and noble metal clusters has fo- 
cused considerable attention along the last years due to the 
different properties detected respect to their corresponding 
bulk-phases. Additionally, these cluster properties have de- 
monstrated to be strongly dependent on the cluster size. On 
the other hand, noble metal clusters, and silver clusters in par- 
ticular, present important applications from the scientific and 
technological point of view in interesting researching fields 
such as catalysis H]|2]|3][4]|5]|6), and optics 1171 151 l9l [TOll . 

The structural and electronic properties of silver clus- 
ters im[l2[32]m[l5j[T6) have been widely studied, however 
unlike the optical properties irTTl IT7l [T81 [20l |2T1 [22l [23l. 
The optical absorption spectra of silver clusters have been 
studied for more than two decades ago. Fedrigo et al. |fT9ll20l 
found that Agjy (N<40) clusters show a simple dominant 
peak in the low energy region, exceptions were attributed to 
the coexistence of low energy isomers. These results were 
confirmed theoretically for small-sized silver clusters (n=5- 
8) llli On the other hand, Idrobo et al. [18| have recently 
carried out calculations of the absorption spectra of Agn for 
different structural isomers, however no full agreement with 
experiments was found. 

On the other hand, experimental optical spectroscopy, aided 
by recent advances in the corresponding theoretical tools, has 
proved to be a powerful means for obtaining information 
about the geometrical and electronic structure of molecules 
and small clusters. Since the structure of a cluster is gene- 
rally not directly accessible experimentally, its characteriza- 
tion must be done with a combination of experiment and the- 
ory. In previous works we have shown how a comparison of 
the measured optical spectrum with theoretical calculations 
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for different isomeric forms of a given cluster can provide 
a useful diagnosis of the geometrical structure of the clus- 
ter l24l [25l . However, conventional optical absorption spec- 
troscopy experiments for studying transition and noble metal 
clusters are difficult due to several technical issues, such as the 
extremely low density of the produced clusters, and theoreti- 
cal spectroscopy is a perfect framework to complement most 
of experimental analysis. 

In this paper, we present calculations of the optical ab- 
sorption spectra of the three most stable isomeric forms of 
the Agn cluster up to 8 eV by using an extremely accurate 
technique for treating the electronic excitations: the time- 
dependent density functional theory (TDDFT) ll26l l27l |281 
under the formalism developed by Casida et al. I29ll30l . 

Our motivations for this work have been the following: 
(a) first, to show a comparison between the present first- 
principles TDDFT results and the available experimental mea- 
surements, in order to improve the description of previous 
theoretical studies by reducing the numerical uncertainty; (b) 
since significant differences have been found among the spec- 
tra of the different isomers of Agn (see Fig.[T|, to propose the 
comparison of experimental and TDDFT absorption spectra 
as a powerful tool to elucidate between the different isomeric 
forms of small clusters; and (c) finally, to predict Agn pho- 
toabsorption spectrum for energies higher than 5.5 eV, a range 
of energies which has not been considered in previous experi- 
mental or theoretical studies for this system (our calculations 
extend the range of excitation energies up to 8 eV). 

II. METHOD 

Prior to the calculation of the photoabsorption spectra, the 
lowest energy geometrical structure of each isomer was deter- 
mined by using the Density Functional Theory ||3D . The ion- 
electron interaction is modelled by replacing the ionic cores 
(ls 2 2s 2 p 6 3s 2 p 6 d 10 4s 2 p 6 d 10 5s 1 Krypton-like core) by pseu- 
dopotentials within the relativistic scheme of Hartwigsen and 
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FIG. 1: View of the three most stable structures of Agn used for the 
calculation of the photoabsorption spectra. The ground-state struc- 
ture is labelled as G-S, and the two closest in energy isomers as ISOl 
and IS02, respectively. 



coworkers [32], previously validated for several metal ele- 
ments in pure and organometallic clusters 041 l25l [33l . This 
pseudopotential is built in such a way that is separable by 
construction, is highly accurate, and has an analytical form 
that can be specified by a small number of parameters. For 
exchange and correlation effects we have used the generalized 
gradient approximation (GGA) 113411 . 

As initial geometries for the structural optimization of the 
Agn isomers were considered the three most stable struc- 
tures previously obtained by Fernandez et al. |[T2l . Starting 
from those geometries we performed structural relaxations 
by making use of the Broyden algorithm P5l with a con- 
vergence criterion of 1CP 2 eV/A in the net forces on every 
atom. The structural relaxations, performed with the ABINIT 
package l36l . did not produce significant variations with res- 
pect to the original structures. In order to test the influence of 
using a well-tailed exchange-correlation potential in the struc- 
tural properties we have additionally checked the simple GGA 
by comparison with a functional with improved asymptotic 
behavior: the van Leeuwen-Baerends (LB94) GGA functio- 
nal 11371 . which led to practically identical optimized struc- 
tures. The geometries obtained, which are then used in the cal- 
culations of the photoabsorption spectra, are shown in Fig.[T] 

Once the ground-state structures of the different isomers 
were established, we performed the calculation of the exci- 
tation spectra. For this purpose, we have used the TDDFT 
formalism [26 27l [28], implemented in the real space code 
octopus [38 1. The theoretical foundations underlying the 
TDDFT calculations, as well as the computational scheme, 
have been presented elsewhere If29l l38l [39l l40l . and here 
we only summarize the main points. The photoabsorption 
spectrum is calculated by using the formalism developed by 
Casida ||29l [30). Application of this technique to obtain 
the oscillator strengths requires a previous calculation of the 



ground state electronic structure, that is, the occupied elec- 
tronic states, and also the unoccupied states. The whole set 
of occupied states, joint to twenty extra unoccupied states for 
each spin channel, were necessary to reach convergence of the 
excitation spectra for an energy range up to 8 eV. After that, 
each excitation peak is broadened by a Lorentzian profile to 
give the photoabsorption cross section 

A 2 

^abs(e)=>J (e _ e . ) 2 +A 2. (D 

where e is the energy, e; are the discrete excitation energies 
obtained by the Casida method, and the value of the parameter 
A determines the full width at half maximum. 

We have performed the calculations of the photoabsorp- 
tion spectra with the GGA, using the Perdew-Burke-Ernzerhof 
parameterization BTI of electronic correlation. The corres- 
ponding occupied and unoccupied single-particle states are 
used in the Casida method with the adiabatic LDA exchange- 
correlation kernel. Although not entirely consistent, this is a 
common practice because the use of a better kernel does not 
significantly affect the results as long as the single-particle or- 
bitals have been calculated with a suitable static v xc (r) poten- 
tial 11251 . This computational framework has been used suc- 
cessfully in the calculation of the optical spectra of atoms and 
clusters |24][25][33][42]|43[4l[45]g^ case 
of free atoms, the energies for one-particle low-energy excita- 
tions agree with experiment to within 5-10 per cent fl31 . For 
clusters of s-p metals, the errors in the position of the absorp- 
tion peaks are usually in the 0. 1-0.2 eV range [45 46 , 47 1, and 
the results are also very reasonable for carbon and BN clus- 
ters Il48"ll4"9ll . and transition metal, organometallic, and molec- 
ular and van der Waals clusters ll24l l25l [331 142|, even in the 
high-energy range. 

This formalism, used in combination with both adiabatic 
local density (ALDA) and GGA approximations for the XC 
effects, has provided successful results, establishing this 
technique for the calculation of excitations in this kind of 
nanoscale systems as an effective tool in the characterization 
of a great number of cluster properties and behaviors, such as 
the elucidation of ground-state structures, the characterization 
of dimeric growths in molecular clusters, and the theoretical 
justification of collective electronic behaviors [24 25, 33, 42 1. 



m. RESULTS AND DISCUSSION 

Fig. [T] shows the calculated equilibrium structures of the 
three lowest-energy isomers of the Agn system, starting from 
the geometries obtained by Fernandez et al. Ifl2l . The results 
do not reveal significant variations respect the original struc- 
tures; where, during the relaxation process, all the atoms in 
the clusters are free to move. The obtained G-S configura- 
tion is based in a seven-atom pentagonal bypiramid, as pre- 
dicted in previous references lfT3"l [181 l22l . The ordering of 
the most stable isomers presented here do not fully agree with 
these previous studies, but is in complete agreement with the 
accurate structural analysis given by Fernandez et al. (for a 
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FIG. 2: Calculated photoabsorption cross sections (in arbitrary units) for the three most stable isomers of Agn (see Fig.|TJ for energies up to 
8 eV. All the absorption strengths are shown in the same scale, and superposed to the photoabsorption profiles vertical dashed lines have been 
included for a better visualization of the different regions of the spectra. 



more detailed description about the structures see Ref. [ 12]). 
However, this factor can be due to the small differences in the 
total energy detected between the different, almost energy- 
degenerated, isomers shown in the recent literature (in the 
range of meV). In the present study the gap in the total energy 
between the G-S structure and the two most stable isomers is 
0.08 and 0.15 eV for ISOl and IS02 structures, respectively. 
Nevertheless, and according to the previous theoretical stud- 
ies, the ground-state of this kind of systems has demonstrated 
to be highly sensitive to the common DFT numerical imple- 
mentations. 

Figure [2] shows the photoabsorption cross sections for the 
three most stable isomers of Agn up to 8 eV as obtained by 
TDDFT under the Casida formalism. The finite width of the 
absorption peaks in an experiment - linked to the accessible 
resolution - is mostly determined by the temperature. In our 
calculations, on the other hand, the width of the peaks is an 
artifact and may be reduced or enlarged by modifying the A 
parameter. We have broadened the peaks in the way shown in 
the figure (with A=0.05 eV), so as to simulate the approximate 
resolution yielded by most of the experiments of photoabsorp- 



tion and photodissociative spectroscopy |24|. 

All three spectra are rich in features in the whole energy re- 
gion up to 8 eV from the begining of the absorption. Three re- 
gions can be clearly distinguished in all the spectra (separated 
by dashed lines superposed to each profile in the picture). A 
region of no absorption at all covers the energy range up to 
2 eV (common to all isomers), a region of strong absorption 
is observed at energies between 2 and 6 eV, and another re- 
gion of weaker absorption at energies larger than 6 eV, with 
an absorption strength dependent on the isomer in particular. 
To our knowledge, no experimental or theoretical study has 
been carried out for Ag n in a range of excitation energies up 
to 8 eV. In the region between 2 and 6 eV, where the absorp- 
tion strength is more pronounced, the spectra are dominated 
by a sharped profile whose position rather varies between the 
different clusters, and formed by pronounced peaks at 3.6, 4.2 
and 4.6 eV. For the ISOl and IS02 isomers an additional peak 
at 5.2 eV is found, not present for the G-S structure. Although 
this wide absorption region is always located between 2 and 6 
eV, and despite of other previous TDDFT theoretical photoab- 
sorption studies in metallic clusters l33l . it cannot be inter- 
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FIG. 3: Calculated optical absorption spectrum of the G-S structure 
of Agn (upper pannel), and three experimental spectra of Agn mea- 
sured in a Ne matrix at 6 K [23], and Ar matrices at 10 [20| and 28 
K 1211 (lower pannels). The computed spectrum has been broadened 
using Lorentzian profiles with a width of 0.05 eV. 



preted as a collective plasmon since it is possible to differen- 
ciate each single-particle excitation forming the profile. It is 
noticeable that the most pronounced peaks forming this strong 
absorption region between 2 and 6 eV are given for the IS02 
isomer, due to the degeneracy of intermediate electronic states 
arising from its high symmetry. In fact, the spectra of the 
other two isomers, where the electronic degeneracy of the in- 
termediate occupied states is much lower, do not present such 
distinct absorption peaks. On the other hand, in the high exci- 
tation energy region (6-8 eV) all the isomers show a first zone 
of moderate absorption up to 7 eV, which is not dependent on 
the particular geometry, and a final part where the strength of 
the absorption increases. As in previous papers 11251 l33ll . we 
propose that at this energy range the high excitation energies 
can be interpreted for this system in terms of the contributions 
between transitions involving strongly hybridated s — d elec- 
tronic orbitals. It is noticeable that this effect is less important 
for the IS02 structure, where, according to its high symmetry 
and Kohn-Sham levels diagram, the sd hybridation is rather 
pronounced. 

The low energy response of these isomers up to 5.5 eV 
can be compared with the experimental results of the Refs. 
1 20 1, 1211 . and 11231 for the measured spectra of medium-size 
silver clusters trapped in noble gas matrices at several tempe- 
ratures. Figure[3]shows the experimental absorption cross sec- 
tions of Agn in a Ne matrix at 6 K, and Ar matrices at 10 and 
28 K (lower pannels), and the TDDFT spectrum for the ob- 
tained Agn G-S structure (upper pannel). The estimated un- 



certainties of the relative cross sections in this kind of experi- 
ments use to be around 10% (the uncertainties in the absolute 
cross sections are even larger). The experimental absorption 
cross sections show a clear three peaks structure, where the 
relative strength of the peaks depend strongly on the tempe- 
rature and the matrix surrounding the trapped cluster, but not 
their excitation energies. The peaks arise at 3.2, 3.7 and 4.4 
eV, with a region of no absorption up to 5.5 eV. The TDDFT 
spectrum, for the G-S isomer, shows a similar three peaks pro- 
file located at 3.6, 4.2 and 4.6 eV. As we can appreciate in 
Fig. [3] there exist a very good agreement between the calcu- 
lated and the experimental peak positions if we just allow for 
a constant shift of around 0.2-0.4 eV. This shift to the red of 
the experimental spectra can be justified by the additional ef- 
fect of the nobel gas matrices to the photoabsorption response, 
which has been quantified by Fedrigo et al. [20], and Sieber 
et al. [ 50 1 in 0.2-0.25 eV. However, most of authors coincide 
in that value of the shift may not be appropriate because ad- 
ditional contributions could arise, such as compression effects 
of the cluster by the matrix, or matrix-induced polarization of 
the cluster charge density [ 1 8 1 . On the other hand, previous 
authors suggest that the deviation in the experimental pho- 
toabsorption spectra of the small silver clusters showing sev- 
eral isolated peaks, instead of only one pronounced excitation 
detected in large-sized clusters lfl9l 1201 . is due to the pres- 
ence of several isomers in the experimental beams |[T8l |2"T1 . 
However, the present study reveal theoretical spectra for three 
different Agn isomers where such isolated peaks are found. 
From these results, we can conclude that the presence of these 
distinct peaks in experimental spectra is not only due to the 
existence of several isomeric forms in the experimental sam- 
ples since it seems to be an intrinsic optical property of this 
kind of aggregates. Other discrepancy sources in the spec- 
trum may be due to additional factors: (i) the error introduced 
in the calculation by the necessarily approximated exchange 
and correlation potential, (ii) the expectable influence of the 
silver clusters coupling to the surrounding noble metal matri- 
ces, and (iii) the possible presence in the experimental sample 
of a mixture of various isomers, coinciding with other authors. 

The calculated spectrum for this isomer provides the best 
matching with experiment, and this confirms the assignment 
of the ground state structure as predicted even by the total 
energies. We then conclude that the calculated spectrum of 
the most stable isomer obtained provides a very accurate de- 
scription of the experimental cross section and permits to con- 
jecture with a high confidence that the isomer present most 
abundantly in the Agn cluster beams has that structure. This 
is reasonable, since this is indeed the calculated lowest energy 
isomer. 



IV. SUMMARY AND CONCLUSIONS 

We have performed calculations of the optical absorption 
cross sections of three structural isomers of Agn. For that 
purpose we have first employed DFT to optimize the geome- 
tries of those isomers. For each of those structures TDDFT 
has been used to calculate the photoabsorption spectrum. In 
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the calculations, we used the generalized gradient approxi- 
mation for the exchange and correlation potential, relativistic 
norm-conserving pseudopotentials, and the Casida formalism 
to calculate the excitation spectrum. This approach yields 
spectra for the different isomers that are significantly diffe- 
rent. The spectrum corresponding to the most stable structure 
obtained matches well the experimental results, and we con- 
clude that this is the isomer predominantly present in the labo- 
ratory cluster beams; confirmation comes from the fact that 
this is the isomer with the lowest calculated energy. The effi- 
ciency of DFT and TDDFT for the calculation of both ground 
state and spectroscopic features of molecules and clusters per- 
mits to use these theoretical methods, as we have shown here, 
as a complement to the experimental analysis, which in many 



occasions lacks some important information, such as the pre- 
cise geometrical atomic arrangement of the probed systems. 
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